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Abstract 

A homogeneous color magnetic field is known to be unstable for the fluctuations 
perpendicular to the field in the color space (the Nielsen-Olesen instability). We 
argue that these unstable modes, exponentially growing, generate an azimuthal 
magnetic field with the original field being in the z-direction, which causes the 
Nielsen-Olesen instability for another type of fluctuations. The growth rate of the 
latter unstable mode increases with the momentum pz and can become larger than 
the former's growth rate which decreases with increasing pz- These features may 
explain the interplay between the primary and secondary instabilities observed in 
the real-time simulation of a non-expanding glasma, i.e., stochastically generated 
anisotropic Yang-Mills fields without expansion. 
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1 Introduction 



Apparent success of hydrodynamic approaches in describing bulk properties 
of the data in relativistic heavy-ion colhsions at BNL-RHIC [1] has stimulated 
great interest in understanding the mechanism for (possible) short-time ther- 
malization of the system from the initial colliding stage. Preceding studies on 
this subject [2] include the scenario based on the plasma instability (in par- 
ticular, the Weibel instability) caused by coupling between hard particles and 
soft fields [3], the 'bottom-up' scenario based on the perturbative scatterings of 
hard particles [4], and so on. In addition to them, we have recently proposed 
a novel scenario based on the Nielsen- Olesen (N-0) instability [5] which is 
characteristic of the configuration of a uniform magnetic field in non-Abelian 
gauge theories [6,7]. This scenario applies the earliest time evolution of the 
system in the heavy ion collision when the separation between hard and soft 
degrees of freedom is not clear. 

In the high-energy limit of collisions, the incident nuclei can be treated as 
two sheets of classical transverse gluonic fields in the Color Glass Conden- 
sate (CGC) framework [8]. This is the QCD effective theory for dense small-x 
degrees of freedom with treating the large-x components as random colored 
sources moving on the light cone. With the initial condition given by CGC, 
it has been shown that longitudinal color electric and magnetic fields are pro- 
duced just after two infinitesimally thin nuclei pass through each other [9]. 
Since the correlation length of the fields on the transverse plane is typically 
~ with Qs being the saturation scale of colliding nuclei, the emerging 

configuration after the collision is made of many longitudinal color flux tubes. 
The produced gluonic system evolving from this unique configuration towards 
locally thermalized state, providing an equilibrated plasma will be formed, is 
called glasma [10]. 

The system of this classical field remains longitudinally boost-invariant dur- 
ing time evolution even though it rapidly expands in the transverse plane and 
is stretched out in the beam direction. Especially, positive longitudinal pres- 
sure is never generated [10,11], and therefore no thermalization is expected to 
be achieved. In reality, however, boost-invariance of the system is naturally 
violated due to finite thickness of the colliding nuclei and quantum fluctua- 
tions [12] both of which are in general rapidity dependent. These rapidity- 
dependent effects may be treated as small fluctuations on top of the boost- 
invariant classical configurations. The first investigation on the consequences 
of rapidity-dependent fluctuations was attempted in a numerical simulation 
of the classical SU(2) Yang-Mills field in an expanding geometry (i.e., in the 
proper-time and rapidity t-t] coodinates) [11] and it was shown that the soft 
fluctuations indeed grow exponentially, contributing positively to the longi- 
tudinal pressure, although the adopted initial fluctuations seemed somewhat 
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too small. This simulation strongly suggests that the boost-invariant system 
is unstable with respect to the rapidity-dependent soft fluctuations. 

Later, analytic studies of the instabilities in the glasma were presented in 
our works [6, 7] within simplified configurations for the background fields. We 
revealed that the characteristic features of the rapidity-dependent fluctuations 
found numerically in Ref. [11] are qualitatively reproduced as the instability 
of the homogeneous color magnetic field a la Nielsen-Olesen [5]. For example, 
the maximum longitudinal momentum of the unstable modes contributing to 
the pressure is shown to increase linearly with the proper-time r in accord 
with the result in Ref. [11]. 

Recently, another unstable behavior in the classical SU(2) gauge field sim- 
ulation was reported in Ref. [13] (see also Ref. [14] for SU(3) case). In this 
simulation, it was observed that (i) there are primary and secondary insta- 
bilities, the former of which is identified with the Weibel instability by the 
authors, and that (ii) the secondary growth rate is larger than the primary 
growth rate and even increases peculiarly with p^. This secondary instability 
was discussed as nonlinear effects of the primary instability in terms of the 
re-summed self-energy diagrams in Ref. [13]. In fact, this simulation is quite 
relevant for physics of the glasma in heavy-ion collisions. The simulation was 
performed in the Cartesian coordinates without expansion (unlike the glasma), 
but adopted an extremely anisotropic initial condition in the momentum space 
which was strongly inspired by the CGC initial condition. More explicitly, the 
initial configuration was stochastically generated according to the distribution 
(in the temporal gauge) 

with the color index a = 1, 2, 3 and the Lorentz index j = x, y, z. The normal- 
ization parameter C is fixed b y the en ergy density. The variance for the 
transverse momentum p_L = y^p^ + may be identified with the saturation 

scale while A-^ for the longitudinal momentum p^ is taken as small enough 
Az ^ A to give a 5(p^)-like distribution practically. We emphasize here that 
the initial configuration thus generated will have a nontrivial correlation over 
the transverse distance ~ l/Qs while it will be very smooth in the z direction^ 
which is quite similar to the longitudinal flux tubes with rapidity-dependent 
soft fluctuations in the glasma [11,7]. Therefore, what was done in Ref. [13] 
essentially corresponds to the numerical simulation of a non- expanding glasma. 

However, there are two points on which we have to be careful in comparing this 
numerical simulation with the expanding glasma. First, the actual simulation 
starts with the supplementary condition = —A'^ = consistently with the 
Gauss law. This means that the initial configuration is purely color magnetic in 
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contrast with the expanding glasma which has initially both color electric and 
magnetic flux tubes. Second, although the magnetic flelds generated by the 
distribution (1) are almost homogeneous in the longitudinal direction, there 
is no preferred direction of the flelds in general. This is in contrast with the 
expanding glasma whose dominant components are longitudinal. Still, it is 
certainly possible for some of the flux tubes to have longitudinal polarization 
in the non-expanding glasma. 

In the present paper, we shall discuss a possibility that the primary and sec- 
ondary instabilities found in Ref. [13] are both induced by the N-0 type in- 
stability studied in [6,7]. The 'primary' instability occurs in a single flux tube 
in which there is a strong magnetic fleld in the longitudinal direction. This 
analysis is essentially equivalent to the case with the expanding glasma [6,7]. 
We will further argue that the 'primary' N-0 instability will induce a color 
electric current in the z direction parallel to the original magnetic fleld, which 
then generates an azimuthal magnetic fleld around the current owing to the 
Ampere law. This new magnetic fleld can become strong enough to cause the 
'secondary' N-0 instability. The "Pz dependence of this subsequent instability 
turns out to be consistent with the flndings in Ref. [13]. 

The present paper is organized as follows. In the next section, we flrst explain 
the N-0 instability of a homogeneous magnetic fleld directed to the z direc- 
tion as an approximation to the magnetic flux tube. We point out qualitative 
similarities with the primary instability found numerically in Ref. [13]. Then 
we argue how the color electric current is induced along the z direction by the 
unstable modes. Stability analysis of the azimuthal color magnetic fleld gen- 
erated by the induced current is presented in Section 3. Section 4 is devoted 
to summary and discussions. 



2 The N-O instability and induced current 



In this section, we flrst review the N-0 instability of a homogeneous magnetic 
fleld in the SU(2) Yang-Mills theory which corresponds to a simplifled situ- 
ation of the magnetic flux tube in the non-expanding glasma. For technical 
reasons, we formulate everything in the Abelian decomposed representation 
as in the original presentation by Nielsen and Olesen [5]. We then discuss the 
consequences of this N-0 instability. 
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2.1 The N-0 instability in the Ahelian decomposed representation 



We study the SU(2) Yang-Mills theory in the temporal gauge Aq = 0. It has 
been well known that the configuration of a homogeneous magnetic field in 
non-Abelian gauge theories exhibits instability with respect to fiuctuations 
perpendicular to the magnetic field in the color space. This is called the 
Nielsen- Olesen (N-0) instability [5]. This can be most easily understood when 
we recast the theory in terms of the U(l) "electromagnetic field" = A^ 
and the "charged vector fields" 0^ = {Aj^ + iA'j^)/\/2 : 

£ _ _ 1 pa paiMU 

= - \Uur - \\D,<Pu - D,<P,\' + igrct>lct>. + - Ml? , (2) 

where we introduced the U(l) field strength f^^ = d^Ay — d^A^ and the 
corresponding covariant derivative = d^j^ + igA^. Notice that we can always 
gauge transform a homogeneous magnetic field to the 3rd color direction so 
that it can be represented by the U(l) gauge field. Then, we treat the U(l) 
gauge field A^ as a background field, while the charged vector fields 0^ as 
small fluctuation^^ Note that the fiuctuations can in principle vary in the z 
direction, while the background field does not. 

We consider a uniform, and time-independent color magnetic field -B > 
in the z direction which can be generated by A^ = {—^By, 0). The 0- 
dependent part of the energy density is now written as (in the temporal gauge 
Ao = h = 0) 

= |0+|2 + + 11)0+1^ + \D(l)-\^ + 2^5|0+|2 - 2gB\(j)-\^ 

+ + - \g\r4>' - 4>'rf (3) 

with Z) = V — igA and 0^ = :^(0^ ± ^0^) in the 'spin' basis. In the second 
line we neglected the surface term and exploited, for simplicity!^ the condition 
Di0* = 0, which is consistent with the non-Abelian part of the Gauss law of 
the SU(2) gauge theory 

+ %gA'<\>' = , (4) 

since the background field is now time-independent A* = 0. Furthermore, the 
physical solution to the Yang-Mills equation must satisfy the Abelian part 
of the Gauss law 

d,)c - ig{^^r - r¥) = , (5) 

which simplifies to cffcff* — = when A^ = 0. We stress here that the 
SU(2) gauge field theory inevitably involves the non-minimal spin coupling 



^ We ignore the fluctuation of the U(l) gauge field since it is stable. 

^ In fact, we are able to obtain the same results without imposing this condition. 



5 



[Df^, D''](j)*^(j)y ~ igf^'^(j)*^(f)i, between and 0^ in addition to the minimal ones 
— ^\D^(f)^\'^ when decomposed with respect to the U(l) subgroup. While this 
non- minimal coupling entails mixing of 0^ and 0^, it can be diagonalized by 
the introduction of 0^ as seen in the second line of Eq. (3). Thus, in a magnetic 
field, the charged vector fields 0^ acquire the "Zeeman" energy ±2gB through 
this non-minimal coupling. This is the origin of the N-0 instability as we will 
see below. 

Within the linear approximation with respect to the fluctuations 0^, we find 
the equation of motion for 0^ 

{d^-D^±2gB)^^ = 0, (6) 

with = — \g'^B'^{x'^ + +'i'gB{ydx — xdy). We can solve it in a similar 
way as in the quantum Hall effects and obtain the eigenfrequencies with the 
longitudinal momentum Pz as 

oj"^ =pI + gB{2n+\m\-m + l)±2gB , (7) 

where = 0, 1, 2, ■ ■ ■ is the principal quantum number of a two-dimensional 
harmonic oscillator and m = 0, ±1, ±2, ■ ■ ■ is the eigenvalue of two-dimensional 
angular momentum = i{xdy —ydx). The factor in the brackets can be writ- 
ten as 2n + \m\ —771 + 1 = 2N + 1 with = 0, 1, 2, ■ ■ ■ specifying the Landau 
levels. Each Landau level is degenerate in m > 0, but we simply set m = 
since it is irrelevant to the discussion below. Note that the frequency for the 
lowest Landau level = n = of 0~ (and 0^* = (0^* + 20^*) /a/2, too) allows 
pure imaginary values for small pz satisfying p^ < gB due to the non-minimal 
coupling, which indicates that the fluctuation grows exponentially 0~ cx: e'''* 
with the growth rate 7 = |Im cul. This is the instability Nielsen and Olesen 
pointed out some time ago [5]. In this ideal case of the uniform magnetic field, 
7 is given by 

7(pJ = |Im ^1 = ^gB-pl. (8) 
Clearly, the fastest instability occurs at = with the maximum growth 
rate 7max = \/gB. Here, we should notice that existence of the unstable mode 
at P2 = is a unique feature of the N-0 instability in clear contrast with 
the Weibel instability, which is familiar in the plasma-instability scenario. 
In fact, the Weibel instability takes place when charged particles with an 
anisotropic distribution move in a soft inhomogeneous magnetic field with 
the 'minimal' particle-field coupling [3]. Inhomogeneity of the soft magnetic 
field is indispensable for an anisotropic distribution of the charged particles to 
show filamentation along this very inhomogeneity inducing a current, which 
then creates a magnetic field additively to the original field. Thus, the Weibel 
instability cannot occur in the homogeneous limit — 0. 

The simplest situation discussed above still has some relevance to the inhomo- 
geneous background magnetic field on the transverse {xy) plane as far as the 
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typical modulating distance R± and the typical strength B of the magnetic 
field satisfy the condition R± ^1/ y/gB. This is because the unstable mode in 
a uniform magnetic field B is localized in a transverse region of a size Xj^fgB 
as is evident from the explicit form of the solutioilEl (in the representation): 

^unstable(^.a:x.P^) 0^ exp ^{j>z)t - ^QBx^}^ . (9) 

As discussed in Introduction, the expanding glasma in heavy-ion collisions is 
initially an assembly of many longitudinal color flux tubes whose transverse 
correlation length is typically given by Q^^ . Since the color magnetic field 
inside the flux tube is roughly estimated as -B ~ 0{Ql/ g), the transverse size 
of the unstable mode of the N-0 instability is also given by \j ^ gB ~ Q~s^ ■ 
Therefore, in this case, the condition > 1/ ^gB is merginally satisfied, 
indicating the relevance of the N-0 instability in the expanding glasma [7]. 
The same should be true for the glasma without expansion if the similar 
estimate for the color magnetic field strength holds. In fact, as mentioned in 
Introduction, the numerical simulation in Ref. [13] was performed choosing the 
initial condition similar to that of the expanding glasma. Therefore, it seems 
plausible that the primary instability found in Ref. [13] is due to the Nielsen- 
Olesen instability. The observed growth rate of the primary instability indeed 
shows very similar dependence as Eq. (8), especially it remains nonzero at 
P. = 0. 

A more rigorous analysis would require to solve Eq. (6) numerically with the 
magnetic field B generated randomly on the transverse plane via Eq. (1). Then 
a negative eigenvalue tu^ of a solution indicates an unstable mode. Although 
existence of an unstable mode in such a generic background is not obvious 
at all, it is quite reasonable that the (largest) growth rate of the unstable 
mode, if exists, will be smaller than the estimate \J~g\B^^\ where |-Bmax| is 
the maximum value of the inhomogeneous magnetic field. This can be intu- 
itively understood if one recognizes that the magnetic field B corresponds to 
a 'potential' in Eq. (6): a positive (negative) B works as an attractive (repul- 
sive) potential well for the fluctuation field 0~ and the unstable modes are 
regarded as 'bound states' trapped in this potential. If one considers a bound 
state (for small pz) in a single potential well with a certain depth, the 'bind- 
ing energy' becomes smaller in general as the potential gets more localized. 
In other words, if we define an effective magnetic field B^q with the (largest) 
negative eigenvalue of u"^ as u"^ = —gBcs, then we have B^s < |-Bmax|- Further- 
more, we may expect that the qualitative dependence of the growth rate on 
the momentum pz will be unchanged from the uniform field case. This is be- 
cause the momentum pz always enters the eigenfrequency u'^ in the same way 
as in Eq. (7), as far as there is no pz dependence in the background magnetic 
field. After all, the pz dependence of the growth rate 7 in an inhomogeneous 

^ This solution satisfies the Abelian part of the Gauss law (5). 
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magnetic field will be represented [6] with the effective magnetic field as 

l{p,)^^gB,^-pl, (10) 
and the maximum momentum for the instability will be given by p^^^ = 



2.2 Induced current and the azimuthal magnetic field 

The exponential growth of the unstable modes cannot last forever. When 
the amplitude of the unstable fluctuation 0~ becomes sizeable, we need to 
include the effects of nonlinear interaction between fluctuations which were 
ignored in the stability analysis of the homogeneous magnetic field. In fact, 
the charged matter fields have a double-well potential similar to the Higgs 
field (see Eq. (3)), and thus it is quite reasonable that the instability will 
cease when the fluctuation becomes as large as 0~ ~ ^JSB/g, corresponding 
to the bottom of the potential. However, this is not the end of the story. In 
this subsection, we discuss the consequences of this enhanced fluctuations. 

We claim that the enhanced fluctuations will induce a large color electric cur- 
rent in the longitudinal direction, which then generates an azimuthal magnetic 
field according to the Ampere law0 Here we define the 'induced U(l) current' 
for the U(l) gauge field in the Lagrangian: = 5C/5A^. Its ^-component is 
given by 

r =Y^ig {(p"*D'(P" - (D"0°)*0"} . (11) 

a=± 

Being made of the charged fluctuations 0^, this current flows in the presence 
of a color electric field in the z direction. This also implies that the current will 
be enhanced if the fluctuations grow exponentially. Recall that the unstable 
mode can have longitudinal momentum ^ \/^^, and thus one may count 
the order of the covariant derivative in the current as ~ 0{^gB). If the 
fluctuation grows up to the same order as its saturation value 0" ~ 0{\J B / g)., 
magnitude of the induced current can be roughly evaluated as 

r^O{g-^- B/g) = 0{{gBf"/g) ~ 0{Ql/g) , (12) 

where we have used an estimate ^JgB ~ Qg. Note that this current can be 
parametrically very large. 



^ Of course, the x and y components of the current are also enhanced by the N-0 
instability, but they generate a magnetic field which helps to cancel the original 
magnetic field in the z direction. 
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The induced current in the z direction creates azimuthal magnetic field B 
around it according to the Ampere law: J = rot B (in the present case 
= -dr{rB^) if one considers only the z component), where the azimuthal 
magnetic field can be as strong as the initial magnetic field: 

B' = 0{Ql/g) . (13) 

For a localized current > in a finite region of the size 1/Qs on the 
transverse plane, the Ampere law implies that the azimuthal magnetic field 
becomes maximum at the brim of the current with the strength estimated as 
in Eq. (13), and it decays only slowly a.s B^{r) oc 1/r outside the region of the 
nonzero current. In the next section, we argue that thus generated azimuthal 
magnetic field is again unstable against the Nielsen-Olesen mechanism. 

Admittedly, our discussion on the induced current (12) presented here is only 
order-of-magnitude estimate. More accurate evaluation of the magnitude will 
require an analysis of nontrivial time evolution of the background field (in 
particular, the longitudinal color electric field E"" = —dtA^) due to coupling 
between the background field and the enhanced fiuctuation, which certainly 
goes beyond the linear analysis. Also more inputs from non-perturbative nu- 
merical simulations will be very useful. One may be able to argue some possible 
mechanisms to generate the induced current, but at present it is quite difficult 
to provide a solid picture for that. We thus leave the definite analysis as an 
open problem for future study. 



3 Effects of azimuthal magnetic field B 

We have seen in the previous section that the homogeneous color magnetic 
field directed to the 3rd color and to the z spatial direction undergoes the 
N-0 instability, and that the enhanced fiuctuations induce a current in the 
z direction, which in turn generates an azimuthal magnetic field B^. In this 
section, we shall investigate the consequences of the presence of this newly 
created magnetic field. As already discussed, the induced current starts to 
fiow according to the generation of longitudinal color electric field E^. Thus, 
the azimuthal magnetic field will become stronger with increasing time. In 
this section, however, we treat B^ as time-independent since we are interested 
in instabilities which will take place in a very short period. Still, it would be 
very helpful to investigate two different situations: (i) when B^ is weak enough 
and can be treated as perturbation to the original magnetic field B^, and (ii) 
when there is a strong B^ . In the latter case, we will ignore the effects of B^ 
which will be screened by the 'primary' N-0 instability. 

Stability analysis in the presence of the azimuthal magnetic field is conve- 
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niently formulated in the cylindrical coordinates (r, 9, z) with x = r cos 9 
and y = r sin 9. For simplicity, we consider a constant color magnetic field 
B = {B'' = Q,B\B') which is generateci3 by the gauge field : {A\A\A') = 
(0, \rB^, —rB^). We again treat the charged vector field 0^ as fiuctuation and 
perform the linear approximation in the equations of motion for 0*. Then one 
finds three coupled equations among cff , (p^ and (f)^ (in the temporal gauge 



AO 



0): 



dl<j)'' - D^(j)'' + -Dg(j)^ + + 2igB'(f)^ - 2igB^(/)' = , (14) 



2 „ \ 



- D^<P' - -De(p' + ^<P" - 2igB'(tf = , (15) 
9,20" - i:>20" + 2i(75V = , (16) 



where = Dl + + \,dr + and Dq = ^de -igA^. Note that it is the 
azimuthal field B^ which couples the field 0" to the transverse field 0*". In fact, 
when the azimuthal magnetic field is absent B^ = 0, Eqs. (14) and (15) reduce 
to Eq. (6) for 0± = {(j)'' ± ^0^)6^^772. 



3.1 B as perturbation to B^ 



Before we present the stability analysis of the azimuthal magnetic field, let us 
examine how the results obtained in the previous section will change in the 
presence of small perturbation of B^. First of all, we rewrite Eqs. (14) and 
(15) in terms of 0^ and 0" (in the Pz representation) as 

[d^ + {p. - gB'rf - \.{rdr) - ± 2gB''^ 0± =2igB'^' . (17) 

Compared with Eq. (6), one indeed verifies that the azimuthal field B^ adds 
the coupling with 0" and shifts the longitudinal momentum in our gauge. 
Consider the unstable fiuctuation 0~ and ignore the stable one 0". When the 
azimuthal magnetic field is much weaker than the longitudinal magnetic field 
B^ B^, we take only the term linearly dependent on B^: 

{{d^ - - 2gB^) - 2gB\pz} 0" = , (18) 

where we have introduced -D^p^ = V^^ ~ pI ~ ^o'^iB^^r'^ for zero angular 
momentum states, with the subscript (0) indicating that B^ = 0. Taking 
the last term in the curly brackets as perturbation, we find that the lowest 
eigenfrequency (squared) is modified in the leading order of perturbation as 

ul^,=pl-gB^-2gB'{r)pz, (19) 
5 = IdgA' - dzA\ = dzA'- - drA^ and = ^drirA'^) - ^deA^ . 
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Fig. 1. A schematic picture of the growth rate 7 = |ImLi;| of the 'primary' N-0 
instabihty. SoUd (dashed) hne is the growth rate with (without) a weak azimuthal 
magnetic field. 



where the average (r) is taken over the unperturbed solution given by Eq. (9), 
i.e., (r) = rdr(p~*r(j)~ / rdr(f)~*(f)~ oc Xj^fgW. This eigenvalue <jj^=q can 
become negative for the momentum such that 

- + gB\r) <p,< + gB'{r) , (20) 

yielding instability with the growth rate 

liPz) - ^gB- + 2gBe{r)p,-pl , (21) 

where we have ignored the contribution quadratic in B^ . The maximum mo- 
mentum for the instability is given by p™^^ ^ ^/qB^ + gB^{r) = y^gB^{l + 
cB^/B^) with c = y^7r/2, and the unstable region in the pz space slightly 
shifts to larger p^. On the other hand, the maximum value of the growth 
rate 7 = |Ima;| is still given by \fgB^ in the leading order perturbation. A 
schematic picture of the growth rate as a function of is shown in Fig. 1. 

Let us summarize here our understanding of the 'primary' N-0 instability. The 
simplest picture with a constant magnetic field B^ should be modified at two 
points: (i) inhomogeneity of the magnetic field on the transverse plane, and 
(ii) the presence of B^ . The first point works to diminish the growth rate, and 
the second one works to shift the instability region in the p^ space slightly to 
the right, to form an oval shape. These two effects seem to be qualitatively in 
agreement with the primary instability obtained in the numerical simulation 
[13]. 
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3.2 Instability of the azimuthal magnetic field B' 



Let us finally present the stability analysis of the azimuthal magnetic field B^ . 
In order to understand the effects of B^ , we consider the simplest case with 
B^ = which is however reasonable because the original magnetic field will be 
weakened by the enhanced fluctuations. Then, Eqs. (14) and (16) are nearly 
diagonalized by combining the fields as (/?^ = (0^ ± z0'')/v^ : 

- D'^+ + ^ (V'^ - V^) + 2^?i?V+ = , (22) 

- i^V" - A (^^ - - 2^?i?V" = , (23) 



where we have again restricted ourselves to the states with dg = and the 
operator D is now given by = D^+d^+^dr. When B^ > 0, the structure of 
these equations suggests that instabilities will occur for ip~ due to the negative 
potential term —2gB^. Therefore, we solve the equation for ip~ assuming that 
ip~^ is small compared to ip~ . By setting ip~^ = 0, one finds an equation for 

I --l-^-J- + - sB'rf + - 2gB'] ^"(r) = ^"(r) . (24) 
V r dr dr ^ ' Ir-^ J 

This is the equation to be solved. Note that the equation has only one pa- 



rameter po = Pz/\JgB^^ which becomes evident if one rewrites it in terms of 



a dimensionless variable p = J gB^r: 



1 d d 2 1 

p dp dp 2p 



P— + (P - Po)' + ^-2}^- =e^- , (25) 



where e = u'^/gB^ is the normalized eigenfrequency. If one further rescales 
the fluctuation field as ^~{p) = p~^^'^ il){p), then one obtains a simple one 
dimensional Schrodinger equation: 

^ +np)V(p)=e^(p), np) = -^ + (p-po)^-2. (26) 



V dp' 4p2 

The explicit form of the potential V{p) is shown in Fig. 2 for three different 
values of po. The potential is singular at p = because of the 'centrifugal' 
potential acting on the m = state, inherent to the vector field, and takes a 
negative value at the minimum which moves outwards with increasing pg. 

Let us first consider the case with large po (high p^). In this case, the centrifu- 
gal potential is negligible around the minimum, and thus the equation reduces 
to that of a harmonic oscillator located at p = po. The eigenvalue of the har- 
monic oscillator then gives approximate solutions to the equation for large po. 
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Fig. 2. Potential V{p) for three different values of pQ. The potential minimum moves 
outwards with increasing po- 



namely, e„ = (2n + 1) — 2 with n = 0, 1, 2, ■ ■ ■ . Clearly, the ground state has 
the negative eigenvalue, e„=o = ~1 because of the last term "—2" whose origin 
can be traced back to the non-minimal coupling between the magnetic field 
and the fluctuation in Eq. (24). The negative eigenvalue u"^ = —gB^ implies 

the existence of an unstable mode with the growth rate 7^ = \fgB^- Notice 
that this result has no Pz dependence. This is because the longitudinal mo- 
mentum Pz (or po) only determines the position of the harmonic oscillator, but 
does not enter the spectrum. The first correction to this harmonic oscillator 
approximation comes from the centrifugal potential, and yields the following 
Pz dependent growth rate (po = Pz/ y qB^) 



7.(P.)^ V^^l-^- (27) 

Therefore, it is an increasing function of Pz, and approaches the harmonic 
oscillator result 7^ = \J gB^ in the limit — * 00. It should be emphasized 
that this instability occurs at larger p^, in contrast to the case of the N-0 
instability in a homogeneous B^ . 

Next we consider the case with small po (low Pz). Notice that the growth rate 
(27) becomes smaller for smaller po. In fact, the equation (26) in the limit 
Po = can be exactly solved by the confluent hypergeometric function with 
the eigenvalues e„ = 4n -|- \/2 for n = 0, 1, 2, ■ ■ ■ [15]. All the eigenvalues are 
positive, and therefore the system is stable when pz = 0. This implies that 
there is a minimum value of Pz for the instability. 
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Fig. 3. Growth rate ^e{Pz) = |Iniw| of the unstable mode in the azimuthal magnetic 
field . Solid line: numerical result from Eq. (26), dashed line: approximate result, 
Eq. (27) valid at large pz- 



Summarizing these analytic studies, we find that the instabihty indeed occurs 
at relatively large Pz with the growth rate given as an increasing function 
of Pzi while there will be no unstable modes for small values of Pz- We can 
directly check that these expectations are indeed the case by solving numer- 
ically Eq. (24). The result is shown in Fig. 3 where the (normalized) growth 

rate 7e/ \l gB^ is plotted as a function of pz/\fgB^. In accordance with the 
analytic result (shown as the dashed line in Fig. 3), the instability appears 

only for larger pz/ \J gB^ > 0.755 and the growth rate indeed increases towards 

the asymptotic value 7^ = \J gB^ with increasing pz- This is the notable fea- 
ture and consistent with the behavior observed for the secondary instability 
in Ref. [13]. We also draw the wavefunctions ip{p) in Fig. 4 corresponding to 
the three different potentials shown in Fig. 2. The wavefunctions are well lo- 
calized around p = po and vanishes at p = 0. The same is true for the original 
fluctuation field <^^(p) = p~^^'^'ip{p). Therefore, we expect that, as far as Pz 
is not too large, our analysis with the homogeneous azimuthal magnetic field 
is not a bad approximation to more realistic cases with B^{p = 0) = and 
S^(p> 1) ~ l/p. 
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4 Summary and discussions 



We have studied instabilities in the configuration with the azimuthal magnetic 
field in the SU(2) Yang-Mills theory, as a possible mechanism for the secondary 
instability observed in the classical stochastic simulation [13]. 

The instability of an expanding glasma numerically found in Ref. [11] was pre- 
viously interpreted as the Nielsen-Olesen instability of the longitudinally uni- 
form color magnetic fields produced in between two sheets of the Color Glass 
Condensates [6,7]. Although the simulation in Ref. [13] was performed in the 
Cartesian coordinates, the adopted initial condition was extremely anisotropic 
in the momentum space, which is inherent to the collision dynamics, and we 
consider the primary instability observed there has the same origin as the one 
in Ref. [11]. 

We have argued that the charged fields grow up due to the N-0 instability 
and are accelerated in the evolving background gauge field, inducing a current 
possibly along the z axis. This current will generates the azimuthal magnetic 
field, whose strength can grow exponentially in time at early stage. The 
induction of the current and its strength is the most speculative part of our 
discussion beyond the linear analysis. 

We have shown that, once the azimuthal magnetic field becomes strong enough. 
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this configuration is subsequently accompanied by the N-0 instability with re- 
spect to another type of fluctuation. We have analyzed this new instability in 
some detail, and shown that the growth rate increases with the momentum 
Pz. This is in contrast with the growth rate of the primary N-0 instabihty 
which decreases with increasing pz and vanishes at some value of pz- The ef- 
fects of the azimuthal magnetic flelds on the primary instability are to enlarge 
the value of the momentum at which the growth rates vanish. These flndings 
coincide quite well with the results in the simulation in Ref. [13]. 

The N-0 instability scenario for the secondary instability is generic to the 
Yang-Mills systems with an anisotropic conflguration, and therefore it should, 
in principle, apply to the case of expanding glasmas, too. In the numerical 
simulation in Ref. [11], however, the secondary instability was not clearly 
observed. This might be because they started the simulation with extremely 
small fluctuations. On the other hand, in Ref. [13], the initial amplitudes of 
the fluctuations were not assumed small (actually, this scale non-separation 
was one of their motivations for taking the numerical approach). Nevertheless, 
it is intriguing to note that, reexamining the pressure evolution started with 
relatively larger fluctuations shown in Ref. [11], one can barely see a subtle 
kink, which may hint the secondary instability. In fact, the increase of the 
longitudinal pressure at late time can be attributed to the generation of the 
transverse magnetic flelds, and this is consistent with our scenario. Another 
indication of the secondary instability is an abrupt increase of the maximum 
longitudinal momentum I'raaxi'r) of the unstable modes. Since the slope C of 
the linear growth t'max(T) = Ct at earlier time is found to be proportional to 
the growth rate [6,7], it is plausible to relate the sudden increase of i^max(T) to 
the emergence of another instability with a larger growth rate. Other possible 
scenarios are previously discussed in Ref. [11]. It is thus a pressing subject 
to analyze speciflcally the course of the expanding glasma evolution from the 
viewpoint of the N-0 instability. 

Our analysis seems to indicate that in the glasma evolution the sequential 
appearance of instabilities play a key role to make the system more random 
and turbulent towards thermalization. Of course, our discussion is based on 
the simple stability analysis and speculation on the induced current. Especially 
more solid discussions on the generation of the azimuthal magnetic fleld and 
the back-reaction from the well-developed unstable mode will require careful 
analyses on the system conflguration and evolution which will be accessible in 
numerical simulations. Moreover, simulations in Ref. [13] should involve the 
Weibel instability mechanism also in principle, although they were formulated 
in terms of the classical flelds. Hence, coexistence of the Weibel and the N-0 
instabilities in the evolving system is worthwhile to be studied. In this context, 
more joint efforts in analytic and numerical approaches are to be accomplished 
in this fleld. 
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